Superintegrability of the Caged Anisotropic Oscillator 

N. W. Evanj] and P. E. VerrieiQ 

Institute of Astronomy, Madingley Rd, 
University of Cambridge, CBS OHA, UK 
(Dated: August 15, 2008) 

Abstract 

We study the Caged Anisotropic Harmonic Oscillator, which is a new example of a superinte- 
grable, or accidentally degenerate Hamiltonian. The potential is that of the harmonic oscillator 
with rational frequency ratio (I : m : n), but additionally with barrier terms describing repulsive 
forces from the principal planes. This confines the classical motion to a sector bounded by the 
principal planes, or cage. In 3 degrees, there are five isolating integrals of motion, ensuring that 
all bound trajectories are closed and strictly periodic. Three of the integrals are quadratic in the 
momenta, the remaining two are polynomials of order 2(1 + m — 1) and 2(1 + n — 1) . In the quan- 
tum problem, the eigenstates are multiply degenerate, exhibiting l 2 m 2 n 2 copies of the fundamental 
pattern of the symmetry group SU(3). 
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I. INTRODUCTION 



The subject of accidental degeneracy has fascinated physicists since the early days of 
the quantum theory. It has long been known that the trajectory in an r-squared force (the 
Coulomb or Kepler problem) is a conic section, so that every bound orbit is an ellipse and 
therefore closed and strictly periodic. The advent of the quantum theory saw the deduction 



of the eigenfunctions for the Hydrogen atom anc 



thence the realization that the bound 



states of the Coulomb problem are degenerate This accidental degeneracy is a 

consequence of a hidden symmetry group SO (4) that is not manifest to the eye. It causes 
all bound trajectories to be closed in the classical problem. 

There is a fundamental connection between accidental degeneracy and the separability 
of the Schrodinger or the Hamilton- Jacobi equations in more than one coordinate system 
and therefore the existence of additional conserved quantities or integrals of motion. This is 
mentioned in a number of the famous texts of the old quantum theory - such as Born's Me- 
chanics of the Atom and Sommerfeld's Atomic Structure and Spectral Lines. For example, 
the Coulomb problem is separable in both spherical polar and rotational parabolic coordi- 
nates. The former leads to the conservation of the angular momentum vector, the latter 
to the conservation of the Laplace-Runge-Lenz vector. The quantum mechanical operators 

n 

close to form the algebra of SO (4) (see, for example, [4| for a review). The accidental de- 
generacy is a consequence of additional integrals of motion Q| and so such systems are often 
called superintegrable |6|]. 

Systematic investigations of all the possible combinations of coordinate systems for which 
the Schrodinger and Hamilton- Jacobi equation can separate have now been carried out j3, [si, 
9[ . In three degrees of freedom, this yields 13 distinct superintegrable systems, all of which 
have classical integrals of motion or quantum operators that are quadratic in the canonical 
momenta and all of which exhibit accidental degeneracy. This includes familiar systems such 
as the Coulomb problem and the isotropic harmonic oscillator. 

Nonetheless, this does not provide a comprehensive explanation of the phenomenon of su- 
perintegrability. For example, in three degrees of freedom, the Hamiltonian of the anisotropic 
harmonic oscillator with rational frequency ratio is 

H = \(PI+P 2 y +Pl) + k(l 2 x 2 + m 2 y 2 + n 2 z 2 ), (1) 
where l,m and n are integers and k is a constant. The Hamilton- Jacobi or Schrodinger 



equations clearly separate in rectangular Cartesians. If I : m : n = 2 : 1: 1, then the 
Hamiltonian also separates in the rotational parabolic and elliptic cylindrical coordinate 



systems giving rise to additional conserved quantities and accidental degeneracy. How 



ever, if l + m+n > 4, then the Hamiltonian is still superintegrable 
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121 ]. even though it 



now only separates in rectangular Cartesians. Further examples of systems which are super- 
integrable but not separable in more than one coordinate system include the Calogero-Moser 
problem {13I . Q] and the generalized Coulomb problem 15 ] . 

The purpose of this paper is to introduce another superintegrable Hamiltonian, closely 
related to ([I]), namely 

H= l -(pl+p 2 y + pt) + k{l 2 x 2 + mY + nV) + % + % + % (2) 

We shall refer to this as "the Caged Anisotropic Oscillator" , as the presence of the barrier 
terms confines the motion to an octant defined by the principal planes. When I — m — n — 1, 
;his becomes the Smorodinsky-Winternitz system, on which there is an extensive literature 



Like the Smorodinsky-Winternitz system, the Caged Anisotropic Oscillator has five in- 
tegrals of motion and it exhibits accidental degeneracy in quantum mechanics. However, 
unlike the Smorodinsky-Winternitz system, the Hamilton- Jacobi and Schrodinger equations 
only separate in rectangular Cartesians. In this paper, we demonstrate that the Hamilto- 
nian (T5]) is superintegrable using the method of projection in §2. Then we discuss both the 
classical and quantum problems in some detail in §3 and §4 respectively. 

II. PROOF OF SUPERINTEGRABILITY 

Let us start with the observation that the commensurate anisotropic oscillator in N 
dimensions possesses 2N — 1 functionally independent integrals of motion, equal to the N 



energies of each individual oscillator and the N — 1 phases differences between them 11]. 
In six dimensions, we have 
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where the rij are positive integers and the Sj are Cartesian coordinates. If we now introduce 
coordinates (x, y, z, 9 X , 9 y , 6 Z ) according to 

s\ = xcos0 x , S2 = xsmO x 
S3 = y cos 8 y , S4 — y sin 9 y 
s 5 = zcos9 z , s 6 = zsin9 z 

and let the frequencies n\ = n 2 = I, n% = = m, = n§ = n, we have the Hamiltonian 



H 



2 \ y x A y z z z 



; /;2 2 1 2 2 1 2 2\ 

k(l x + m y + n z ) 



(4) 



Since the coordinates (0 x ,8 y ,8 z ) are ignorable we can set the conjugate momenta to con- 
stants. Making the substitutions p$ x = ki, Pg y = k 2) p 2 dz = k 3 gives us the Caged Harmonic 
Oscillator Hamiltonian (TSl). 



For the Hamiltonian in 6 dimensions given by ([3]), every bound trajectory is closed. 
Similarly, in the reduced 3 degrees of freedom Hamiltonian (j2J), every bound trajectory is 
also closed and the system is still superintegrable. 



k, = 0.1 k 2 =0.2 k 3 =0.3 k = 0.5 






FIG. 1: A typical orbit for a frequency ratio of I : m : n = 1 : 1 : 3. 
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FIG. 2: A typical orbit for a frequency ratio ofZ:m:ra = l:2:3. 



III. CLASSICAL MECHANICS 
A. The Integrals of Motion 

The Hamiltonian ([2]) clearly separates in rectangular Cartesian coordinates to give the 
first three integrals of motion as the three energies of oscillation 



1 „2 , /_2„,2 , k 2 



I x = L - v \ + kl 2 x 2 + ^ (5) 



h = -P 2 y + km 2 y 2 + -1 (6) 
2 y 



As the system inseparable in these coordinates, it is easy to see that the trajectories of the 
orbits are (c.f. 



x2 = 4k + (m? ~ m) <-(^«-<»» («) 
^ = 4k + (iJp " J*) V2 cos( ^ m e " ( °> + Ci > < 9) 
22 = Ji + (4^-^) 1/!cos ' VSn ( t - t »» +C2 » < 10 » 

where to and the c, are constants. The remaining two integrals are the phase differences 
between the orbits, say c\ and c^- If we say too that \m — l\ < \n — l\ < In — ml, then the 



integrals are also of lowest order possible in the momenta. First, let us define 

2 

£ = — — = cos(v / 8A?/(t - t )) (11) 

v 2 - t— 

rj = n = cos(V8km(t - 1 ) + c x ) (12) 

2 

C = = cos{V8kn{t - £ ) + c 2 ) (13) 



where a = h/(2l 2 k), (3 = I 2 /(2m 2 k) and 7 = I 3 /(2n 2 k) and 
A = 



Al A k 2 l 2 k) ' \Am 4 k 2 m 2 k 

C= ML_M 1/2 (14) 

The derivation of both integrals is similar and will be demonstrated with the case of cj. The 
first phase difference is given by (c.f., the discussion of the anisotropic oscillator in |12| ) 

m 

c\ = arccos?] — - arccos £ (15) 

Taking the cosine gives 

cos(/ci) = cos(/ arccos rj) cos(m arccos£) + sin(Z arccos rj) sin(m arccos£) 

= T^TUO + ^^ 2 T;(V)TL(0 (16) 



where 7] and T m are the Chebyshev polynomials of the first kind 18] and T( and T' m are 
their derivatives with respect to the arguments rj and £. The time derivatives £ and fj are 
equal to 2xp x /A and 2yp y /B respectively. It is more convenient to express the integral as 

h = {2k) l+rn A m B l cos(/ Cl ) (17) 

which is of order 2(1 + m) in the momenta, but can be reduced to order 2(1 + m — 1) since 
the two highest powers of the momenta can be removed though a combination of the energy 
integrals. 

The corresponding integral for the second phase difference c 2 is 

J 5 = (2k) l+n A n C l cos(Zc 2 ) (18) 

where 

cos(/c 2 ) = 7](C)T n (0 + gJ^/(Cra£) (19) 



which can be reduced to order 2(1 + n — 1) in the momenta. It is easy to verify that both 
I4 and I5 are integrals of motion by showing that the Poisson bracket with the Hamiltonian 
vanishes. They are also functionally independent, as may be verified by computing the rank 
of the appropriate Jacobian. 




FIG. 3: A typical orbit for a frequency ratio of 2 : 3 : 4. 
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FIG. 4: A typical orbit for a frequency ratio of 2 : 3 : 4 with no potential barriers (k\ = k<i = k% = 
0). The dotted lines show the field of view of the corresponding orbit with the potential barriers 
as shown in Fig [3l 

B. The Group Theoretic Approach 

Rodriguez et al. [19( have also recently examined this system, and derived the classical 
integrals. Ingeniously, they look for invariants under SO (2) x 50(2) x SO (2) corresponding to 
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the transformation (J4]). They then search for combinations of these invariants that commute 
with the Hamiltonian ([3]) under the Poisson bracket. Such quantities will also necessarily be 
integrals of the reduced Hamiltonian of the Caged Anisotropic Oscillator. Rodriguez et al 
find integrals of motion that are rational functions in the momenta, but here we show how 
to adapt their method to give integrals that are polynomial in the momenta. 
Let us start by introducing the complex variables 

Z\ — pi — iiV2ksi, zi = P2 — iiV2ks 2 , 
Z3 = Ps — imV2ks3, -24 = Pa — imV2ks4, 

Zb = Pb — in\/2ks5, z§ = p§ — inV2ksQ, (20) 
so that the Hamiltonian ([3]) is just 

H = \j^\z,?. (21) 

i=l 

Now, following Rodriguez et al, we look for invariants under the generators of rotations in 
the (31,^2), (53,54) and (55,^) planes. For the (si,S2) plane, they include 

zl + zj, zl + 4, \ Zl \ 2 + \z 2 \ 2 . (22) 

with similar results holding for the (53,54) and (55, sq) planes. Expressions like \zi\ 2 + 
\z2 1 2 clearly commute with the Hamiltonian ([3]) and are just the separable energies in the 
oscillation in the coordinate directions. The remaining quantities do not commute with ([3]), 
but it is possible look for an invariant that is a function of the two expressions z\ + z\ and 
z\ + z\ that does. Therefore, we require that 

{H 6 J(z 2 1 +z 2 2 ,z 3 2 + z A 2 } = 0. (23) 

Inserting H 5 from ([3]), this gives the complex invariant 

R = (% + 4) m (4 + 4Y (24) 

whose real part 

R + R = {zl + z 2 2 ) m {zl + ztY + {z\ + zlT{zl + zlY (25) 



s 



is a polynomial of order 2(£ + m). Modulo an unimportant overall numerical factor, it is 
the same as the polynomial invariant found earlier in eq (ITT)) . Similarly, the invariant (ITS)) 
is just 

$ + -lr(^ + 4Y + (4 + 4T(4 + zlY m 

up to a numerical factor. 




k, = 0.1 k 2 = 


0-2 k 3 =0.3 k = 0.5 
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i x = [1,2,3] »=[3,2,1] 
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FIG. 5: A typical orbit for a frequency ratio of 3 : 4 : 5. 




FIG. 6: A typical orbit for a frequency ratio of 3 : 4 : 5 with no potential barriers (ki = ki = = 
0). The dotted lines show the field of view of the corresponding orbit with the potential barriers 
as shown in Fig [5j 



C. The Orbits 



It is interesting to plot out the orbit of a particle in the potential. As the Hamiltonian 
is sup erintegr able, all bound orbits must be closed curves. Using a standard Bulirsch-Stoer 
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integrator [2Cj to solve the equations of motion, some example orbits are plotted. These are 
shown for various frequency ratios in Figs [U [2J [3] and [51 The orbits are confined to a box, 
defined by the limits a±A, and similar. Figs H] and El show the corresponding cases to Figs 
Eland El bu t with no potential barriers. Here, the trajectories are the well-known Lissajous 
figures 2l), and it can be seen how the orbit in the general case is a reflection and slight 
distortion in the x, y and z axes. 



IV. QUANTUM MECHANICS 

The Schrodinger equation is separable in rectangular Cartesians, and reads: 

3 



8=1 



(27) 



where a>j are the integer multipliers of the frequencies, corresponding to /, m, n in the previous 
section, and h — 1. The separable solution has wavefunction ^ = n^LiVv where the 
individual wavefuntions are (c.f., 

I22|) 

where N Ui is the normalisation constant given by 

N nt = w 1 t /2 (2k) 1 ^J(2wfk)^T(n t + 1)/T(m + l±u t ) 



(28) 



(29) 

and L v n . are associated Laguerre polynomials, the T are Gamma functions and Vi = |(1 + 
8fcj) 1//2 . The quantised energy is given by 

/ I 7/. \ 

(30) 



E = 2V2kJ2wi (rn + l±if) 
»=i ^ ' 



The degeneracy of each energy level with quantum number N = W\n\ + w^n^ + ^3^3 is 
therefore the same as that of the three dimensional anisotropic harmonic oscillator with 
rational frequency ratio w\ : W2 '■ W3 (listed for example in [^j) In the simplest case, if the 



frequency ratio is 1 : 1 : n then the degeneracy is given by 

g(N) = 



N 
n 



I IliV + l-^ 



N 
n 



(31) 



where [N/n] denotes the integer part of N/n. The allowed states for three degrees of freedom 
and the frequency ratio 1:1:2 are shown in Fig. [71 
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To look at the group structure, the annihilation and creation operators can be constructed 



as 



h 
b] 





1 




s/2k 




1 


4wr 


s/2k 



q 2ki d 2 

2wiV2kxi- h 2w 2 kx 2 ^ + WiV2k + — j 

OXi ' xf oxf 

—2wiV2kxi— \-2WjKXj ^ — WiV2k + — -k 

OX; ' xf 0X7 



which annihilate and create quanta of energy in the i direction, that is 



[H,bi] = -2wiV2kbi 
[H,b\] = 2wiV2kb\ 

and, representing ip n . as \rii), act in the following way 



bi\rii) = y/m(ni ± Vi)\ni - 1) 

b\\ni) = y/(m + l)(ni±i/i + l)\ni + l) 



(32) 
(33) 



(34) 
(35) 



(36) 
(37) 



which are identical to those for the Smorodinsky-Winternitz system given by 16|] . As such 
the number operator given by n» = | ([£>*, b\] =F z/j — 1) can be used to again to construct the 
operators 

T^l^^^Mn,^)-^} (38) 
which close under commutation 



(39) 



and give the Lie algebra u(3). 

In the case of the isotropic harmonic oscillator, it is well known that the degeneracy of the 
iVth energy level is (N+l)(N+2)/2, which corresponds to the dimensions of the irreducible 
representations of SU (3) . Even though the symmetry group of the anisotropic harmonic 
oscillator is also SU(3), it is no now longer the case that the degeneracy levels follow the 



pattern 1, 3, 6, 10, 15.... This was already noted as a complication by Jauch & Hi 



there have been a number of possible resolutions proposed in the literature 
Following the lines of argument put forward in [24j, we define 



23J, 



24 



1 
25 



Hi and 



H0 



h\ = rii mod (1V2W3), n 2 = n 2 mod (wiWs), n 3 = 113 mod (1V1W2) (40) 
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from which it follows that 



nx — nxw 2 w 3 + n, n 2 = n 2 wxw s + r 2 , n 3 — n 3 WxW 2 + r 3 



(41) 



where < rx < w 2 w^, < r 2 < wxw% and < r 3 < wxw 2 . This divides the energy levels 
into w\w\w\ subsets according to the values of rx,r 2 and r%. From eq. fl30l) . the energy levels 
become 



E = 2V2kJ2 



{wxW 2 W3)ni + Wi{ri + - =b — 



(42) 



so that the energy levels within each subset (rx, r 2 , r 3 ) have the characteristic degeneracy of 
SU(3). This is illustrated by the color coding in Fig. [7J 




n«, 3 



FIG. 7: Energy levels for the frequency ratio 1:1:2. Objects belong to different sets of (rx,r 2 ,r 3 ) 
values are shown in different color. For each color, the degeneracies are the dimensions of the 
irreducible representations of SU(3). 



V. SUMMARY AND CONCLUSIONS 



The Caged Anisotropic Harmonic Oscillator is a new superintegrable Hamiltonian, 
namely 



// = \(pl +P 2 y + V\) + Hl 2 x 2 + mY + nV) + lL+-j± + l±. ( 13) 



x 2 y 2 z 2 
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If the frequency multipliers are integers, then the Hamiltonian is superintegrable. We have 
found the five isolating integrals for the classical motion in three degrees of freedom. Three of 
the integrals of motion - the energies in each oscillation - are quadratic in the canonical mo- 
menta and arise from separation of the Hamilton- Jacobi equation in rectangular Cartesians. 
The other two integrals are still polynomial in the momenta, but now of order 2(7 + m — 1) 
and 2(7 + n — 1) respectively. If / = m = n = 1, the Hamiltonian becomes the well-studied 



Smorodinsky-Winternitz system 



281 ]. and all the integrals are then quadratic and 



arise from separability of the Hamilton- Jacobi equation. 

The system is interesting for at least three reasons. First, from the perspective of integra- 
bility, there are still very few systems known with integrals of motion that are polynomials in 



the momenta of higher order than 2 



29j . Systematic searches for Hamiltonian systems with 



higher order polynomial invariants have been performed, confirming the impression that 



they are rare 30j, |3jj . Given this sketchy and disparate information, we have no unifying 



theory of the conditions for the existence of such integrals of motion 

Second, from the perspective of superintegrability, if the integrals of motion are all 
quadratic in the momenta, then a classification theorem exists and all systems in flat space 
have been found [3, Sis]. Such systems always arise from separability of the Hamilton- Jacobi 
equation in more than one coordinate system. However, the Caged Anisotropic Oscillator 
joins the Toda Lattice and the Generalized Kepler Problem as an example of a system for 
which some of the integrals are cubic polynomials or higher, and then the superintegra- 
bility does not arise from separability in more than one coordinate system. It would be 
interesting to classify such systems and find all examples in flat space. In particular, the 
Caged Anisotropic Oscillator is the second superintegrable Hamiltonian to be deduced by 
the method of projection introduced in [ijj]. Essentially, the idea is to view superintegrable 
motion in three degrees of freedom as a projection of a higher dimensional superintegrable 
system, such as the Coulomb or Kepler problem, or the harmonic oscillator. Are there any 
more such systems to be found? 

Third, from the perspective of group theory in quantum mechanics, the proper interpre- 
tation of the symmetry or degeneracy group remains unclear. Already in 1940, Jauch & 
Hill [lOj noted that the quantum mechanical problem of the anisotropic oscillator presents 
problems which leaves its symmetry group in doubt. Since that day, there have been a num- 
ber of different suggestions in the literature as to the proper interpretation of the symmetry 
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group 



24 



25l . l26l . [27|. Although these procedure seem reasonable, they are more along the 



lines of a posteriori justification than compelling argument. 
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